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We consider one single copy of a mixed state of two qubits and investigate how its entanglement 
changes under local quantum operations and classical communications (LQCC) of the type p' ~ 
(A ® B)p(A ® B)^ . We consider a real matrix parameterization of the set of density matrices 
and show that these LQCC operations correspond to left and right multiplication by a Lorentz 
matrix, followed by normalization. A constructive way of bringing this matrix into a normal form is 
derived. This allows us to calculate explicitly the optimal local filterin operations for concentrating 
entanglement. Furthermore we give a complete characterization of the mixed states that can be 
purified arbitrary close to a Bell state. Finally we obtain a new way of calculating the entanglement 
of formation. 
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Entanglement of two separated quantum systems im- 
plies that there are non-local correlations between them. 
This feature of non-locality has found practical applica- 
tions in quantum information theory (see for example 
0). Most applications require that both parties share 
maximally entangled states. A realistic preparation and 
transmission of entangled states however yields mixed 
states. Therefore Bennett et al. j|] proposed a protocol 
which allows to obtain asymptotically a nonzero num- 
ber of maximally entangled pure states by carrying out 
collective measurements on a large number of copies of 
entangled mixed states. Their scheme however required 
that the fidelity of the mixed states exceeded 1/2. The 
Horodecki's subsequently showed how mixed states of ar- 
bitrary fidelity could be purified by applying first a fil- 
tering operation on each copy separately Linden et 
al. Q then asked the question if it is possible to obtain 
singlets out of mixed states by allowing only local oper- 
ations on each copy separately. While this is possible for 
pure states, they proved that this is impossible in general 
for mixed states UH, as the best state one can obtain 
is a Bell diagonal state |J. The Horodecki's however 
gave an example of a mixed state that could be puri- 
fied arbitrary close to a singlet state by a process called 
quasi-distillation Q. 

We shed new light on those results by observing that 
filtering operations on two qubits correspond to Lorentz 
transformations on a real parameterization of their den- 
sity matrix. Using Lorentz transformations, this real pa- 
rameterization can be brought into one of two types of 
normal forms, thus giving a characterization of all states 
that can be transformed into each other by local opera- 
tions. Our scheme also yields a new way of calculating 
the entanglement of formation j^] , with as a by-product 
a simple proof of the necessity and sufficiency of the par- 
tial transpose criterion of Peres ||[l0| ■ The main result of 
this letter however is the fact that we provide a construc- 
tive way of finding the optimal POVM's for concentrating 



the entanglement. We show that there exist two classes of 
states corresponding to the two normal forms : the ones 
that can be brought into Bell diagonal form leaving the 
rank of the density matrix constant, and the ones that 
can asymptotically be brought into Bell diagonal form 
with lower rank. This last class contains a subclass of 
mixed states that can be purified arbitrary close to the 
singlet state. 

In this letter we will consider the filtering operations 
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(A g B)p(A g B)^ 
Tr((A <g> B)p(A (g) i?)t 



(1) 



where A' A < I2, B^B < I2. As a local projective 
measurement destroys all entanglement, we will only 
consider the cases det(A) ^ and det(£?) 7^ 0. Let 
us now calculate how the entanglement of formation 
(EoF) changes under these local operations. The EoF 
of a two qubit system can be calculated as a convex 
monotonously increasing function of the concurrence Q . 
As shown in JlT| , the concurrence of p is given by 
max(0, ti — T2 — T3 — T4) with {7^} the singular values 
of X T (a y o- y )X with p = XX^. Under the filtering 
operations we have X' = (A® B)X/y/Tr(A^A <g> B^Bp). 
As {A®B) T (o- y ®(T y ){A®B) = det(A)det(B)(a y ®o- y ), 
this proves the following theorem: 

Theorem 1 Under the filtering operations the 
currence changes as 



con- 



C = C 



|det(A)||dct(B)| 
Tr(AtA <g) BtBp) ' 



(2) 



It will turn out very usefull to introduce the real and 
linear parameterization of the density matrix |p_ 2[| 



(3) 



1 



where the summation extends from to 3 and with o~o 
the 2x2 identity matrix and cr±, o~i, 03 the Pauli spin ma- 
trices. Below we will often leave out the normalization 
of p and R. Note that normalization of R is very simple 
since i?o,o = Tr(p). 

We will now prove how R transforms under the LQCC 
operations (Q). 

Theorem 2 The 4x4 matrix R with elements Rij — 
^ Tr(p(<7^ <8>CTj)) transforms, up to normalization, under 
LQCC operations my as 



R 1 — LaRL^q 



(4) 



where La and Lb are proper orthochronous Lorentz 
transformations given by 



L A =T(A® A*)Ty\det(A)\ 
L B =T{B®B*)Ty\det{B)\ 
/ 1 1 \_ 
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This theorem can be proven by introducing the matrix 
Pkl,k'l' — Pkk',W an d noting that R = 4TpT T . It is easy 
to check that under the LQCC operations (|l|) p trans- 
forms as p' = {A <g> A*)p(B (g) B*) T . Therefore R trans- 
forms as R' = L A RL%\det(A)\\det(B)\ with L A = T(A® 
A*)Tt/|det(A)|, L B = T(B(g> B*)Ty\det{B)\. Using the 
identities Aa y A T = det(A)a y and T^MT* = -o y ® a y 
with M the matrix associated with the Lorentz met- 
ric M — diag[l, — 1, — 1, — 1], it is easily checked that 



LaML a — M — LbML b . Furthermore the determi- 



nant of La and Lb is equal to +1, and the (0, 0) element 
of L is positive, which completes the proof. □ 

As the complex 2x2 matrices with determinant one in- 
deed form the spinor representation of the Lorentz group, 
there is a 1 to 2 correspondence between each La and 
Aj dct(A). It is interesting to note that when both A 
and B arc unitary, the theorem reduces to the well-known 
fact p2| that the rows and columns of R transform under 
50(3), which is indeed a subgroup of the Lorentz group. 

With the above theorem in mind, a natural question 
is to find a decomposition of R as R = L\YaL\ with E 
diagonal and L\, L2 proper orthochronous Lorentz trans- 
formations. This would be the analogue of a singular 
value decomposition but now in the Lorentz instead of 
the Euclidean metric. 

Theorem 3 The 4x4 matrix R with elements R^ — 
Tr(p(ii ® tjj) can be decomposed as 

R = LxZLl (8) 

with L±,L2 proper orthochronous Lorentz transfor- 
mations, and E either of diagonal form E = 



diag[s , si, s 2 , S3] with s > 
the form 



si > s 2 > 



(a b 
d 
-d 
\c a+c-b ) 



I S3 1, either of 



(9) 



with a, 6, c, d real. 

The proof of this theorem is quite technical. It heavily 
depends on results on matrix decompositions in spaces 
with indefinite metric p3[ . We first introduce the matrix 
C = MRMR T which is M-selfadjoint. Using theorem 
(5.3) in |13| ], it follows that there exist matrices X and 
J with C = X~ 1 JX, J consisting of a direct sum of real 
Jordan blocks and XMX T — Nj with Nj a direct sum 
of symmetric nxn matrices of the form [Sij] = ±[<5i+j,„+i] 
with n the size of the corresponding Jordan block. Us- 
ing Sylvester's law of inertia, there exists orthogonal Oj 
such that Nj = O^MOj. It is then easy to check that 
OjX — L\ is a Lorentz transformation. Therefore the 
relations C = MRMR T = ML^lOjJO^Lf hold. Mul- 
tiplying left by M, Sylvester's law of inertia implies that 
there exist a matrix E with the same rank as J such 
that MOjJO T j = EME T . Therefore we have the rela- 
tion RMR T = LiEAfE 7 L\ . If R has the same rank as 
RMR T , this relation implies that there exists a Lorentz 
transformation L2 such that R = L\£jL\. 

Let us now investigate the possible forms of E. As 

Nj = T jMOj has signature (H ), J can only be 

a direct sum of the following form: 4 lxl blocks; 1 or- 
thogonal 2x2 block and 2 lxl blocks; 1 2x2 Jordan block 
and 2 lxl blocks; 1 3x3 Jordan block and 1 lxl block. 
Noting the eigenvalues of C as {A^}, it is easy to verify 
that a "square root" E in the four cases is respectively 
given by 

1. E = diag[vW, VM VP^P with P a 

permutation matrix pcrmutating the first column 
with one other column: 



2. E = diag 



3. E = diag 



4. E = diag 



cos(0) sin(</>) 
sin(0) — cos(^) 



a b 
c a + c — b 



( a 

b Va 2 + b 2 

q — ab a 2 



with 



a = y/\X \ and b = -l/y/2\X \. 

Now we go back to the relation R = LJY.L2- Li and 
L 2 can be made proper and orthochronous by absorbing 
factors —1 into the rows and colums of E yielding E'. 
Theorem (2) now implies that this E' corresponds to an 
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unnormalized physical state, which means that p' corre- 
sponding to £' has no negative eigenvalues. It is easy to 
show that this requirement excludes cases 2 and 4 of the 
possible forms of S. The third case corresponds to (||). 
Furthermore in the first case the permutation matrix has 
to be the identity and |Ao| > max(|Ai|, | A2 1 , | A3 1 ) . Mul- 
tiplying left and right by proper orthochronous Lorentz 
transformations, the elements {si} of this diagonal £ can 
always be ordered as So > Si > s% > |«3 j . 

The case where the rank of C is lower then the rank 
of R still has to be considered. This is only possible if 
the rowspace of R has an isotropic subspace Q for which 
QMQ T — 0. Some straightforward calculations reveal 
that the only physical states for which this hold have 
normal form (^|) with a = b = c and d = or a — b and 
c = d = 0. This completes the proof. □ 
The two normal forms can be computed very efficiently 
by calculating the Jordan canonical decomposition of 
C = MRMR T and of C = MR T MR. It is easy in- 
deed to show that for example in the case of diagonaliz- 
able R the eigenvectors of C form a Lorentz matrix, and 
I Si I = ^J\i(C). Note that we always order the diagonal 
elements such that s > s i > s 2 > In- 
states that are diagonal in R correspond to (unnor- 
malized) Bell-diagonal states with ordered eigenvalues 



Ai = (s + si + s 2 - s 3 )/4 
A 2 = (s + si - s 2 + s 3 )/4 
A3 = (so - si + S2 + s 3 )/4 
A 4 = (s - si - s 2 - s 3 )/4, 



(10) 

(11) 
(12) 
(13) 



whereas states of type (§) correspond to the rank defi- 
cient states 



/ a + c d \ 



b-c 

\ d a-b J 



(14) 



For both cases it is easy to calculate the entangle- 
ment of formation analytically, respectively given by |1J] 
C = max(0, (Ai - A 2 - A 3 - A 4 )/(A X + A 2 + A 3 + A 4 )) = 
max(0, (— s + si + s 2 — s 3 )/(2si)) and C — max(0, \d\/a). 

Let us now consider an arbitrary state p with corre- 
sponding R. Combining theorem (1),(2) and (3), it fol- 
lows that the concurrence of p is equal to the concurrence 
of the state corresponding to £ multiplied by i? o- We 
have therefore proven: 

Theorem 4 Given a state p and associated with this 
state R = L\LL\ , then the concurrence of p is given by 
C = max(0, (— s + si + s 2 — s 3 )/2) or by C = max(0, |<i|) 
depending on the normal form £. 

We thus have obtained a new method of calculat- 
ing the entanglement of formation of a system of two 



qubits. Interestingly, it turns out that this character- 
ization relates the concepts of entanglement of forma- 
tion and of partial transposition Let us therefore 
define Rf? — Tr(p PT <Ti o-j), which changes the sign of 
the third column of R. In the case of diagonal normal 
form of R, it is readily verified that the normal form of 
R PT equals that of R except for the last element where 
s PT = — s 3 . Retransforming £ PT to the p PT -picture, we 
see that the corresponding Bell-diagonal partial trans- 
posed state has minimal eigenvalue so — si — s 2 + s 3 . We 
readily recognize the expression of the concurrence of the- 
orem (||) and therefore this eigenvalue is negative if and 
only if the concurrence exceeds 0. Moreover we know 
that p PT is related to this Bell-diagonal state by some 
similarity transformation A®B which cannot change the 
signature of a matrix due to the inertia law of Sylvester. 
In the case of normal form (^), an analogue reasoning 
shows that p PT has a negative eigenvalue if and only if 
\d\ > 0, which again is necessary and sufficient to have 
entanglement. This completes the proof of: 

Theorem 5 Given a system of two qubits, this state is 
separable if and only if its partial transpose has a negative 
eigenvalue. 

Although this result was already proven by Horodecki 
pol l, we believe the previous derivation is of interest as 
it connects the entanglement measures concurrence and 
negativity. Using this formalism, it indeed becomes pos- 
sible to prove that the concurrence always exceeds the 
negativity, and it is furthermore possible to find a com- 
plete characterization of all states with maximal and min- 
imal negativity for given concurrence [|l5|| . This is im- 
portant because in the two qubit case the negativity is a 
measure of the robustness of entanglement against noise. 

Next we want to solve the problem of finding 
the POVM such as to have a non-zero chance to 
produce a new state with highest possible entangle- 
ment. From equation (|^), the maximum EoF is ob- 
tained with A, B minimizing the expression Tr(j4T^4 £g> 
B t Sp)/(|det(A)det(S)|. Absorbing the factors |det(A)| 
and |det(J3) I into A and B, it is sufficient to consider 
A and B with determinant 1. In the i?-picture, the 
optimization is then equal to minimizing the (0,0) ele- 
ment of R = Li£L^ by appropriate La, Lb- Absorb- 
ing L\ and L 2 into L' A = LaL^M and L' B = LbL^M, 
this is equivalent to finding the optimal vectors I a and 
Ib such that I^YjIb is minimized under the constraints 
l\Ml A = 1 = 1%MI B . 

Let us first consider the case of diagonal £ with el- 
ements so > si > s 2 > |s 3 |. Parameterizing l A as 
(-^/l + ||x|| 2 , x) and Ib as (\A + ||y|| 2 , y), the following 
inequalities hold: l^Yl B > s ^l + ||x|| 2 ^l + \\y\\ 2 - 
s i \\x\\ || y| > so- Therefore the concurrence will be max- 
imized for x — y = 0, leaving £ into diagonal form. 
Collecting the previous results, it follows that if R is di- 
agonalizable, the state with maximal concurrence that 
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can be obtained from it by single copy LQCC operations 
is the one corresponding to £ which is a Bell diagonal 
state. This is in complete accordance with the results 
of Kent et al. [|| . The optimal A and B are thus given 
by the 2x2 matrices corresponding to LjM and L^M. 
The optimal POVM can then be obtained by dividing A 
and B by their largest singular value such that A^A < 1 
and B^ B < 1, fo llowed by cal culating t he square roots 
A c = y/ 12 — A^A and B c = yi-j — B^B which are rank 
1. The optimal POVM's to be performed on the two 
qubits are then given by {A, A c } and {B, B c } respec- 
tively. Note that the probability of measuring (A, B) is 
given by the inverse of the gain in concurrence divided 
by the product of the largest singular values of A and B, 
and that the rank of the Bell diagonal state is equal to 
the rank of the original state. 

If £ is of the form (||) however, things get more 
complicated. An analogous reasoning as in the di- 
agonal case leads to the conclusion that I a and Ib 
are vectors associated with the Lorentz transformations 
bringing (||) into diagonal form. This is however only 
possible in the limit where Ia and Ib contain factors 
rim t ^oo[-\/i + i 2 ,0,0,t] and lim^oo [Vl + i 2 , 0, 0, -t] re- 
spectively. This indeed allows to bring R asymp- 
totically into diagonal form with diagonal elements 
given by [•*/ (a — b)(a + c), d, —d, yj(a — b)(a + c)] and 
off-diagonal elements of order 1/t 2 , yielding a state in- 
finitesimally close to a Bell diagonal state. The probabil- 
ity to get this state during a measurement of the optimal 
POVM however scales as lim^oo 1/t 2 . This is equiva- 
lent to the quasi-distillation protocol by Horodecki Jl0| . 
In this limit of t — > 00 the rank of the new state is less 
than the original one, and its concurrence is given by 
\d\/^(a-b)(a + c). 

In the case where a — 6 = a + c = \d\ we are there- 
fore able to create a state arbitrary close to the singlet 
state. Therefore the only mixed states that can be quasi- 
purified to the singlet state by single copy LQCC oper- 
ations are the rank two states having normal form (|^) 
with a — b — a + c = \d\. 

In conclusion, we obtained new insight into the prob- 
lem of local filtering on one copy of two qubits by intro- 
ducing the notion of Lorentz transformations on a real 
matrix parameterization of their density matrix. This 
matrix can be brought into one of two types of normal 
forms. These normal forms contain all the information 
about the entanglement of formation and reveal an ele- 
mentary connection between concurrence and the partial 
transpose criterion of Peres. Moreover, this new formal- 
ism enabled us to derive in a constructive way the optimal 
local filtering operations for concentrating entanglement 
on an arbitrary mixed state of two qubits. This could 
be of great interest in constructing optimal distillation 
protocols. We showed that states of the first type can 



be locally transformed into a Bell diagonal state of the 
same rank with finite probability, whereas states of the 
second kind can asymptotically be transformed into Bell 
diagonal states with lower rank. This last class is of spe- 
cial interest as is contains the mixed states that can be 
transformed arbitrary close to the singlet state. 
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